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ABSTRACT 

^^ I We measure the angular bispectrum of the cosmic microwave background (CMB) 

(^ ' radiation anisotropy from the COBE Differential Microwave Radiometer (DMR) four- 

^N ■ year sky maps. The angular bispectrum is the harmonic transform of the three-point 

^ ■ correlation function, analogous to the angular power spectrum, the harmonic transform 

'"^ . of the two-point correlation function. First, we study statistical properties of the bis- 
pectrum and the normalized bispectrum. We find the latter more useful for statistical 
analysis; the distribution of the normalized bispectrum is very much Gaussian, while the 

'T^ I bare bispectrum distribution is highly non-Gaussian. Then, we measure 466 modes of 

1/^ ' the normalized bispectrum, all independent combinations of three-point configurations 

^^ ■ up to a maximum multipole of 20, the mode corresponding to the DMR beam size. By 

[~^ ■ measuring 10 times as many modes as the sum of previous work, we test Gaussianity of 

^^ . the DMR maps. We compare the data with the simulated Gaussian realizations, finding 

^^ . no significant detection of the normalized bispectrum on the mode-by-mode basis. We 

r"| ! also find that the previously reported detection of the normalized bispectrum is con- 

Mh| sistent with a statistical fluctuation. By fitting a theoretical prediction to the data for 

O ' the primary CMB bispectrum, which is motivated by slow-roll inflation, we put a weak 

-*— * ' constraint on a parameter characterizing non-linearity in inflation. Simultaneously fit- 

Cd ■ ting the foreground bispectra estimated from interstellar dust and synchrotron template 

K^ . maps shows that neither dust nor synchrotron emissions significantly contribute to the 

^ . bispectrum at high Galactic latitude. We conclude that the DMR map is consistent 

H ! with Gaussianity. 

Subject headings: cosmology: observations - cosmic microwave background ~ early 
universe 
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1. INTRODUCTION 

Why study non-Gaussianity of the cosmic microwave background (CMB) radiation anisotropy? 

Inflation (Guth 1981; Sato 1981; Albrecht and Steinhardt 1982; Linde 1982) predicts Gaussian 
primordial perturbations in quantum origin (Guth and Pi 1982; Hawking 1982; Starobinsky 1982; 
Bardeen et al. 1983), implying that two-point statistics such as the angular power spectrum, C;, 
specify all the statistical properties of the CMB anisotropy. Inflation has passed several challenging 
observational tests; the recent CMB experiments (Miller et al. 1999; de Bernadis et al. 2000; Hanany 
et al. 2000) have shown that the universe is flat as predicted by inflation with a fluctuation spectrum 
consistent with an adiabatic scale-invariant fluctuation. 

Several authors have attempted to measure non-Gaussianity in CMB using various statistical 
techniques (e.g., Kogut et al. 1996b); as yet no conclusive detection has been reported except 
for measurement of several modes of the normalized CMB bispectrum on the COBE Differential 
Microwave Radiometer (DMR) sky maps (Ferreira et al. 1998; Magueijo 2000). The existence of 
non-Gaussianity in the DMR data is controversial. If the CMB sky were non-Gaussian, this would 
challenge our simplest inflationary model. 

The angular bispectrum, Bi^i^i^^ is the harmonic transform of the three-point correlation func- 
tion. We carefully distinguish the normalized bispectrum, Bi^i^iJ [Ci^Ci^Ci^) ' , from the bispec- 
trum, Bi^i^i^. Ferreira et al. (1998) have measured 9 equilateral {h = I2 = I3) modes of the normal- 
ized bispectrum, Bi^i^i.^/ [Ci^Ci^Ci^) ' , on the DMR map, claiming detection at l\ = l2 = lz = 16. 
Their result has been under extensive efforts to confirm its significance and origin. Bromley and 
Tegmark (1999) claim that a few individual pixels in the DMR map are responsible for the most 
of the signal. Banday et al. (2000) have proposed an eclipse effect by the Earth against the COBE 
satellite as a possible source of the signal. Magueijo (2000) has measured other 8 inter-/ modes 
of the normalized bispectrum such as Bi^m^i/ (Q-iQCi+i) ' , and claims that scatter of the 
normalized bispectrum among 8 modes is too small to be consistent with Gaussian. Sandvik and 
Magueijo (2000) further report measurement of 24 other inter-Z modes for different lags in /, and 
conclude they are consistent with Gaussian. 

Hence, until now 41 modes of the normalized bispectrum have been measured on the DMR 
map. Here, we simply ask: "how many modes are available in the DMR map for the bispectrum?" 
The answer is 466, up to a maximum multipole of 20 that corresponds to the DMR beam size; thus, 
it is conceivable that the claimed detection of the normalized bispectrum at h = h = h = 16 would 
be explained by a statistical fluctuation, as 9 modes are expected to have statistical significance 
above 98% out of 466 independent modes even if CMB is exactly Gaussian. In this paper, we 
measure 466 modes of the CMB bispectrum on the COBE DMR sky maps, testing the claimed 
detection of the bispectrum and non-Gaussianity. We take into account the covariance between 
these modes due to the Galactic cut, which has not been done in the previous work. 

On the theoretical side, several predictions for the CMB bispectrum exist. Several authors (Falk 



et al. 1993; Luo and Schramm 1993; Gangui et al. 1994) have predicted the primary bispectrum (or 
equivalently three-point correlation function) on the DMR angular scales from slow-roll inflation 
models. Komatsu and Spergel (2001) have extended the prediction down to arcminutes scales using 
the full radiation transfer function. 

In addition to the primary one, secondary sources in the low-redshift universe and foreground 
sources produce the bispectrum through their non-linearity. Luo and Schramm (1993) and Spergel 
and Goldberg (1999) have calculated the secondary bispectrum arising from non- linear evolution 
of gravitational potential; Goldberg and Spergel (1999) and Cooray and Hu (2000) have calculated 
the one from the gravitational lensing effect coupled with various secondary anisotropy sources. 
Komatsu and Spergel (2001) have calculated the foreground bispectrum from extragalactic radio 
and infrared point sources. While the bispectrum is not the best tool for detecting the signature 
of rare highly non-linear events, e.g., textures (Phillips and Kogut 2001), it is sensitive to weakly 
non- linear effects. 

Having theoretical predictions is a great advantage in extracting physical information from 
measurement; one can fit a predicted bispectrum to the data so as to constrain parameters in a 
theory. Since the DMR beam size is large enough to minimize contribution from the secondary and 
the extragalactic foreground sources, the only relevant source would be the primary one. In this 
paper, we fit a theoretical primary bispectrum (Komatsu and Spergel 2001) to the data. 

The Galactic plane contains strong microwave emissions from interstellar sources. The emis- 
sions are highly non-Gaussian, and distributed on fairly large angular scales. Unfortunately, pre- 
dicting the CMB bispectrum from interstellar sources is very difficult; thus, we excise the galactic 
plane from the DMR data. We model the residual foreground bispectrum at high galactic latitude 
using foreground template maps. By simultaneously fitting the foreground bispectrum and the 
primary bispectrum to the DMR data for three different Galactic cuts, we quantify the importance 
of the interstellar emissions in our analysis. 

This paper is organized as follows. In § 2, we define the angular bispectrum, and show how 
to compute it efficiently from observational data. In § 3, we study statistical properties of the 
bispectrum and the normalized bispectrum. We then measure the normalized bispectrum on the 
COBE DMR four- year sky maps (Bennett et al. 1996), testing Gaussianity of the DMR map. In 
§ 4, we fit predicted bispectra to the DMR data, constraining parameters in the predictions. The 
predictions include the primary bispectrum from inflation and the foreground bispectrum from 
interstellar Galactic emissions. Finally, § 5 concludes. In the appendix, we derive the relations 
between the angular power spectrum and bispectrum on the incomplete sky and those on the full 
sky. 



2. ANGULAR BISPECTRUM 

The CMB angular bispectrum consists of a product of three harmonic transforms of the CMB 
temperature field. For Gaussian fields, expectation value of the bispectrum is exactly zero. Given 
statistical isotropy of the universe, the angular averaged bispectrum, Bi-^^i.^.^, is given by 

R — \^ I ^^ ^2 ^^1 (^\ 
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where the matrix denotes the Wigner-3j symbol. The harmonic coefficients, aim, are given by 

aim = [ d'n ^^YCm (n) , (2) 

where r^obs denotes a solid angle of the observed sky. Bi-^i^i.^ satisfies the triangle condition, \li — lj\ < 
Ik ^h+ Ij for all permutations of indices, and parity invariance, ^i + ^2 + ^3 = even. 

We rewrite equation (1) into a more computationally efficient form. Using the identity. 
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we rewrite equation (1) as 

Bhl2h={ 000) J ^ ^hi^)(^l2i^)eh{'^), (4) 

where the integral is not over Jlobs) but over the whole sky; e;(n) already encapsulates the infor- 
mation of incomplete sky coverage through aim- Here, following Spergel and Goldberg (1999), we 
have used the azimuthally averaged harmonic transform of the CMB temperature field, e;(n), 



ei{ri.) = \l^j^^aimYim{n). (5) 



Similarly, we write the angular power spectrum, Ci, as 

^^ = /^«^(-)- (6) 

ei (n) is thus a square-root of Q at a given position of the sky. 

Equation (4) is computationally efficient, as we can calculate e/(n) quickly with the spherical 
harmonic transform for a given I. Since the HEALPix pixels have the equal area (Gorski et al. 
1998), the average over the whole sky, J(i^n/(47r), is done by the sum over all pixels divided by 
the total number of pixels, N~^ Yli ■ 



3. MEASUREMENT OF BISPECTRUM ON THE COBE DMR SKY MAPS 

3.1. The data 

We use the HEALPix-forniatted (Gorski et al. 1998) COBE DMR four-year sky map, which 
contains 12,288 pixels in Galactic coordinate with a pixel size 1?83. We obtain the most sensitive 
sky map to CMB by combining 53 GHz map with 90 GHz map, after coadding the channels A 
and B at each frequency. We do not subtract eclipse season time-ordered data; while Banday et 
al. (2000) ascribe the reported non-Gaussianity to this data, we will argue in this paper that the 
claimed detection of the normalized bispectrum at h = h = h = 16 (Ferreira et al. 1998) can also 
be explained in terms of a statistical fluctuation. 

We reduce interstellar Galactic emissions by using three different Galactic cuts: the 20° cut, 
the extended cut (Banday et al. 1997), and the 25° cut. We then subtract the monopole and the 
dipole from each cut map, minimizing contaminations from these two multipoles to higher order 
multipoles through the mode-mode coupling. The coupling arises from incomplete sky coverage. 
This is very important to do, for the leakage of power from the monopole and the dipole to the 
higher order multipoles is rather big. We use the least-squares fit weighted by the pixel noise 
variance to measure the monopole and the dipole on each cut map. 

We measure the bispectrum, Bi^i^i^^ on the DMR sky maps as follows. First, we measure 
aira using equation (2). Then, we transform aim for —l<m<l into e/(n) through equation (5). 
Finally, we obtain Bi^i^i^ from equation (4), arranging /i, I2, and ^3 in order of h < I2 < hj where 
we set the maximum ^3 to be 20. In total, we have 466 non-zero modes after taking into account 
l^i — lj\ < Ik ^ h + Ij and ^1 + ^2 + ^3 = even. Measurement of 466 modes takes about 1 second of 
CPU time on a Pentium-Ill single processor personal computer. 



3.2. Monte Carlo Simulations 

We use Monte-Carlo simulations to estimate the covariance matrix of the measured bispectrum. 
Our simulation includes (a) a Gaussian random realization of the primary CMB anisotropy field 
drawn from the CO-B-B-normalized ACDM power spectrum, and (b) a Gaussian random realization 
of the instrumental noise drawn from diagonal terms of the COBE DMR noise covariance matrix 
(Lineweaver et al. 1994). For computational efficiency, we do not use off-diagonal terms as they are 
smaller than 1% of the diagonal terms (Lineweaver et al. 1994). 

We generate the input power spectrum, Q, using the CMBFAST code (Seljak and Zaldarriaga 
1996) with cosmological parameters fixed at Jlcdm = 0.25, Oa = 0.7, il^ = 0.05, h = 0.7, and n = 1; 
the CMBFAST code uses the Bunn and White (1997) power-spectrum normalization. 

In each realization, we generate aim from the power spectrum, multiply it by the harmonic- 
transformed DMR beam, Gi (Wright et al. 1994), transform Giaim back to a sky map, and add an 



instrumental noise realization to the map. Finally, we measure 466 modes of the bispectrum from 
each realization. We generate 50,000 realizations for one simulation; processing one realization 
takes about 1 second, so that one simulation takes about 16 hours of CPU time on a Pentium-Ill 
single processor personal computer. 



3.3. Normalized bispectrum 

The input power spectrum determines the variance of the bispectrum. Off-diagonal terms in 
the covariance matrix arise from incomplete sky coverage. When non-Gaussianity is weak, the 
variance is given by (Luo 1994; Heavens 1998; Spergel and Goldberg 1999; Gangui and Martin 
2000) 

(Bli.i,) = {Ci,) {Ci,) {CQ Ai,i,i„ (7) 

where A^^/j/^ takes values 1, 2, or 6 for all Vs are different, two are same, or all are same, respectively. 
The brackets denote the ensemble average. 

The variance is undesirably sensitive to the input power spectrum; even if the input power 
spectrum were slightly different from the true power spectrum on the DMR map, the estimated 
variance from simulations would be significantly wrong, and we would erroneously conclude that 
the DMR map is inconsistent with Gaussian. It is thus not a robust test of Gaussianity to compare 
the measured bispectrum with the Monte-Carlo simulations. 

The normalized bispectrum, Bi-^i^i^/ (Ci^Ci^Ci^) ' , is more sensible quantity than the bare 
bispectrum. Magueijo (1995) shows that the normalized bispectrum is a rotationally invariant 
spectrum independent of the power spectrum, as it factors out fluctuation amplitude in aim, which 

1/2 

is measured by C^ . By construction, the variance of the normalized bispectrum is insensitive to 



the power spectrum, approximately given by Ai^i^i.^. 

One might wonder if the normalized bispectrum is too noisy to be useful, as the power spectrum 
in the denominator is also uncertain; however, we find that the variance is actually slightly smaller 
than Al-^^l^l^. Figure 1 compares the variance of the normalized bispectrum, (^Bf^i^^, / (Ci-^Ci^Ci.^)), 
with that of the bispectrum, {Bf^i,) / {{Ci^) (Ci^) {Ci^)). The top-left panel shows the case of full 
sky coverage. We find that the variance of the normalized bispectrum is precisely 1 when all /'s are 
different, while it is slightly smaller than 2 or 6 when two /'s are same or all Ts are same, respectively. 
This arises due to correlation between the uncertainties in the bispectrum and the power spectrum, 
and this correlation tends to reduce the total variance of the normalized bispectrum. The rest 
of panels show the cases of incomplete sky coverage. While the variance becomes more scattered 
than the case of full sky coverage, the variance of the normalized bispectrum is still systematically 
smaller than that of the bare bispectrum. The normalized bispectrum is thus reasonably sensitive 
to non-Gaussianity, yet it is not sensitive to the overall normalization of power spectrum. 

What distribution does the normalized bispectrum obey for a Gaussian field? First, even for a 



A 



CO 
CVJ 



PQ 
V 




2345671234567 
<B2 , >/(<C, ><C, ><C, >) 



'r2^3 



Fig. 1. — Comparison of the variance of the normahzed bispectrum, \Bf , ,/ (Ci-^Ci^Ci.^)), with 
that of the bare bispectrum, (i^f^^j/- ) / ii^h) i^h) (^^3))- The top-left panel shows the case of 
full sky coverage, while the rest of panels show the cases of incomplete sky coverage. The top- 
right, bottom-left, and bottom-right panels use the 20° cut, the extended cut, and the 25° cut, 
respectively. 



Gaussian field, the probability distribution of a single mode of Bi^i^i^ is non-Gaussian, characterized 
by a large kurtosis. Figure 2 plots the distributions of 9 modes of -Bzijj's drawn from the Monte- 
Carlo simulations (solid lines) in comparison with Gaussian distributions calculated from r.m.s. 
values (dashed lines). We find that the distribution does not fit the Gaussian very well. Then, 
we examine distribution of the normalized bispectrum, Bi^i^i^/ [Ci^Ci^Ci^) ' . We find that the 
distribution is very much Gaussian except for h = h = h = 2. Figure 3 plots the distributions of 
the 9 modes of the normalized bispectrum (solid lines) in comparison with Gaussian distributions 
calculated from r.m.s. values (dashed lines). The distribution fits the Gaussian remarkably well; 
this motivates our using standard statistical methods developed for Gaussian fields to analyze 
the normalized bispectrum. We could not make this simplification if we were analyzing the bare 
bispectrum. Furthermore, the central limit theorem implies that when we combine 466 modes the 
deviation of the distribution from Gaussianity becomes even smaller. 

Ferreira et al. (1998) claim detection of the normalized bispectrum at h = h = h = 16; 
Magueijo (2000) claims that the scatter of the normalized bispectrum for li = I2 — I and 1^ = 12 + 1 
is too small to be consistent with Gaussian. The former has analyzed 9 modes, while the latter has 
analyzed 8 modes. In the next section, we analyze 466 modes, testing the statistical significance 
of the non-Gaussianity with much more samples than the previous work. We calculate Ci from 
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equation (6), and then divide Bi-^i^i^ by {Ci^Ci^Ci.^ ' to obtain the normalized bispectrum. 



3.4. Testing Gaussianity of the DMR map 

We characterize statistical significance of the normalized bispectrum as probability of the mea- 
sured normalized bispectrum being greater than those drawn from the Monte-Carlo simulations. 
We define the probability P as 

p. - ^'" '^'" '^ = / dx FMC(x), (8) 

^V total -^-l^f^'^l 

where /„ is the normalized bispectrum, A'totai = 50, 000 is the total number of simulated realizations, 
and a = 1, 2, 3, 4,..., 466 represent {hMM) = (2,2,2), (2,3,3), (2,2,4), (3,3,4),..., (20,20,20), 
respectively, with satisfying li < I2 < h, Ik — lj\ ^ h ^ k + ^j, and ^i + ^2 + ^3 = even. F^'^{x) 
is the p.d.f of the simulated realizations for the normalized bispectrum, x = I^ ■ The p.d.f is 
normalized to unity: J^^ dxF^ {x) = 1; thus, Pa lies in < Pq, < 1. 



By construction, the distribution of Pa is uniform, if the DMR map is consistent with the 
simulated realizations, i.e., Gaussian. We give the proof as follows. By rewriting equation (8) as 
P^ = /(jl^^^l), we calculate the p.d.f of P„, G(P„), as 

/oo 
dy5[Pa = f{\y\)]F^^\y) 
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Fig. 2. — Distribution of the bispectrum drawn from the Monte-Carlo simulations for the 20° 
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Galactic cut (solid lines). Bi^i^i^/ {{Ci^) (Ci^) {Ci.^)) ' is plotted, where the brackets denote the 
ensemble average over realizations from the Monte-Carlo simulations. The dashed lines plot Gaus- 
sian distributions calculated from r.m.s. values. Each panel represents a certain mode of {h,l2,h) 
as quoted in the panels. 
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Fig. 3. — Distribution of the normalized bispectrum, Bi-^i^i^/ {Ci^Ci^Ci.^ ' , drawn from the Monte- 
Carlo simulations for the 20° Galactic cut. The meaning of the lines is the same as in figure 2. 
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where F^ (y) is the p.d.f of the measured normalized bispectrum on the DMR map, y = I^ . 
Our goal is to see if F^ (y) is consistent with the DMR data being Gaussian. It follows from 
equation (9) that G{Pa) = 1, when F^^^{y) = F^^(y), regardless of the functional form of 
F^^(y). In other words, the distribution of Pa is uniform, if the distribution of the measured 
normalized bispectrum is the same as the simulated realizations. Since our simulation assumes 
the DMR map Gaussian, the P distribution, G{P), tests Gaussianity of the DMR map. If the P 
distribution is not uniform, then we conclude the DMR data to be non-Gaussian. 

A Gaussian field gives equal number of modes in each bin of P. For example, it gives 46.6 
modes in AP = 10% bin: 466xG(P)AP = 466x0.1 = 46.6. If we detect the normalized bispectrum 
significantly, then we find that G{P) is not uniform, but increases rapidly as P increases. 

The top panel of figure 4 plots the P distribution for the three different Galactic cuts. We find 
that the distribution is uniform, and the number of modes in the bin (AP = 10%) is consistent 
with the expectation value for Gaussian fluctuations (46.6). 

To further quantify how well it is uniform, we calculate the Kolmogorov-Smirnov (KS) statistic 
for the P distribution in comparison with the uniform distribution. The bottom panel of figure 4 
plots the cumulative P distribution, for which we calculate the KS statistic. The probability of the 
distribution being uniform is 6.7%, 73%, and 77% for the three Galactic cuts, respectively. 

We have confirmed that the normalized bispectrum at h = h = h = 16 has P = 97.81% for the 
20° cut, P = 99.97% for the extended cut, and P = 99.27% for the 25° cut, as similar to Ferreira 
et al. (1998); however, our result shows that statistical fluctuations explain the significance. We 
conclude that the properties of the normalized bispectrum of the DMR map are consistent with 
CMB being a Gaussian field. 



4. MODEL FITTING 

In this section, we fit predicted CMB bispectra to the measured normalized bispectrum. The 
predictions include the primary bispectrum from infiation and the interstellar foreground bispec- 
trum from the Galactic emissions. Then, we constrain a parameter characterizing the primary 
bispectrum. 
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Fig. 4. — P distribution (Eq.(8)). P is the probability of the CMB normaUzed bispectrum, 
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measured on the COBE DMR 53 + 90 GHz sky map, being larger than 



those drawn from the Monte-Carlo simulations. There are 466 modes in total. The thick dashed, 
solid, and dotted lines represent the three different Galactic cuts as quoted in the figure. The thin 
solid line shows the expectation value for a Gaussian field. The top panel shows the P distribution, 
while the bottom panel shows the cumulative P distribution, for which we calculate the KS statis- 
tic. The KS statistic gives the probability of the distribution being consistent with the expectation 
for Gaussianity as 6.7%, 73%, and 77% for the three Galactic cuts, respectively. 
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4.1. Primary bispectrum 

For the primary bispectrum from inflation, we consider weakly non-Gaussian adiabatic per- 
turbations generated through non-Hnearity in slow-roll inflation. The simplest weak non-linear 
coupling gives 

$(x) = $l(x) + /nl [^l(x) - (^l(x)>] , (10) 

where the square bracket denotes the volume average, and <1'l(x) is a linear Gaussian part of 
curvature perturbations. We call /nl the non-linear coupling parameter, following Komatsu and 
Spergel (2001). 

Salopek and Bond (1990, 1991) and Gangui et al. (1994) show that slow-roll inflation gives 
this coupling; Pyne and Carroll (1996) shows that the second-order general relativistic perturbation 
theory gives this. The former predicts /nl as a certain combination of slope and curvature of a 
inflaton potential (<1>3 = — 2/nl in Gangui et al. (1994); a^ = /nl in Verde et al. (2000)). The 
latter predicts /nl ^ C'(l). Komatsu and Spergel (2001) have given the exact form of Bi-^i^i^ for 
this model; we do not repeat it here, /nl is the parameter that we try to constrain by measuring 
the CMB bispectrum. 

Since the theoretical bispectrum assumes full sky coverage, we must correct it for the bias aris- 
ing from incomplete sky coverage. We use an approximate correction factor for the bias, r2obs/4vr, 
which we have derived in the appendix. Moreover, the theoretical bispectrum must also be con- 
volved with the DMR beam. We use the harmonic transform of the DMR beam, G;, given in Wright 
et al. (1994). Hence, we relate the observed bispectrum to the theoretical bispectrum through 



^M^/3 = ^^MrGz.Gz,G.3- (11) 



Note that r2obs/4vr = 1 — sin|6cut| for an azimuthally symmetric cut within certain latitude 6cut; 
^obs/4'/r = 0.658, 0.577, and 0.5 for |5cut| = 20°, 25°, and 30°, respectively. For the extended cut, 
Oobg/47r = 0.638. 



4.2. Foreground bispectra from interstellar emissions 

Although we cut a fraction of the sky to reduce interstellar emissions from the Galactic plane, 
there should be some residuals at high Galactic latitude. Kogut et al. (1996a) have found signifi- 
cant correlation between COBE DMR maps at high Galactic latitude and COBE Diffuse Infrared 
Background Experiment (DIRBE) map which mainly trace dust emission from the Galactic plane. 

The interstellar emissions are highly non-Gaussian. For example, the one-point p.d.f of the all- 
sky dust template map (Schlegel et al. 1998) is highly skewed. We find the normalized skewness, 
((AT)^) / <^(AT)^) ~ 51. Since these non-Gaussian emissions would confuse the parameter 
estimation of the primary CMB bispectrum, we take the effect into account. 
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We estimate the foreground bispectra from interstellar sources by using two foreground tem- 
plate maps. One is the dust template map of Schlegel et al. (1998); the other is the synchrotron 
map of Haslam et al. (1981). Both maps are in the HEALPix format (Gorski et al. 1998). 

We extrapolate the dust map to 53 GHz and 90 GHz with taking into account spatial vari- 
ations of dust temperatures across the sky (Finkbeiner et al. 1999). We then cross-correlate the 
extrapolated maps with the DMR maps to confirm that the extrapolation is reasonable. We find 
that while the dust-correlated emission in the DMR 90 GHz map is consistent with the extrapo- 
lated dust emission, that in the DMR 53 GHz map is much larger than the extrapolated one. This 
is consistent with the anomalous microwave emission of Kogut et al. (1996a). To take the excess 
emission into account, we multiply our extrapolated 53 GHz maps by factors of 3.66, 2.58, and 2.59 
for the 20° cut, the extended cut, and the 25° cut, respectively. The correction factor is notably 
larger for the 20° cut rather than the extended or 25° cut . This could possibly be attributed 
to the region around Ophiucus which has a free-free spectrum. Note that we do not essentially 
need the correction for the excess emission, as it does not alter spatial distribution of the emission. 
Nevertheless, we do it for convenience of subsequent analyses. 

We also extrapolate the synchrotron map to these two bands, assuming the spectrum of the 
source, T{v) ex v~'^'^. We do not need the extrapolation of the synchrotron template map either, as 
the extrapolation does not alter spatial distribution of the emission in contrast to the dust template 
map in which the extrapolation does alter it. We find no significant correlation between the DMR 
maps and the extrapolated synchrotron maps at both 53 GHz and 90 GHz. 

After coadding the extrapolated 53 and 90 GHz maps with the same weight as used for the 
DMR maps, we measure -B^^/jia from the maps for the three different Galactic cuts, multiplying it 
by Gi^Gi^Gi^ to take into account the DMR beam. 



4.3. Constraints on non-linearity in inflation 

We fit simultaneously the primary, dust, and synchrotron bispectra to the measured bispectrum 
on the DMR map. We use the least-squares method based on a x^ statistic defined by 

x\h) - E f 4°'^^ - E f^^i] (^")..' U^"" - E /.^^') • (12) 

aa' \ i J \ i / 

la is a model bispectrum divided by ((C'/ ) {Ci ) {C^ )) , where j represents a certain com- 
ponent such as the primary, dust, and synchrotron, fj is a fitting parameter for a component j, 
where /primary = /nl is the non-linear coupling parameter (Eq.(lO)). /dust and /sync characterize 
amplitude of the foreground bispectra. 

Caa' is the covariance matrix of the normalized bispectrum, which we calculate from the 
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Monte-Carlo simulations: 



N 



^-' - ]v3T E V^""^'^ - i^^'')) i^^'''^'' - (4^^>) ' (13) 



j=l 



where N = 50, 000 is the total number of realizations. The bracket denotes an average over all 
realizations, (/^^) = N~^ ^^ la ■ Here, we have implicitly assumed the non-Gaussianity weak, 
so that we calculate the covariance matrix from Gaussian realizations. 

As we have observed in the previous section, distribution of la is very much Gaussian; thus, 
X^ifj) should obey the x^ distribution to good accuracy. Hence, minimizing x^ifj) with respect to 
fj gives the maximum-likelihood value of fj as a solution to the normal equation. 



i 

where 



ElUC~')aa'I. 



DMR 
a' 



(14) 



^..-E^"(^"')..'^«'- (15) 

aa' 

We estimate statistical uncertainties of the parameters using the Monte-Carlo simulations; we 
obtain parameter realizations by substituting /^^ for /^^^ in equation (14). 

Figure 5 plots the measured values of the non- linear coupling parameter, /nl, as well as the 
simulated realizations, for the three different Galactic cuts. The measured values are well within 
the cosmic variance: we place 68% confidence limits on /nl as |/nl| < 1-6 x 10^, 1.5 x 10^, and 
1.7 X 10'^, for the 20° cut, the extended cut, and the 25° cut, respectively. 

Figures 6 and 7 plot constraints on /dust and /syno respectively. There is no indication of 
either component contributing to the measured bispectrum significantly. 



4.4. Null test of the normalized bispectrum 

Using x^ defined by equation (12), we can test Gaussianity of the DMR map. While the 
minimization of x^ifj) gives constraints on the parameters, a value of x^ifj) t^lls us goodness-of- 
fit; x^(0) tests a hypothesis of the bispectrum being zero. When x^(0) is either significantly greater 
or smaller than those drawn from the simulations, we conclude that the DMR map is inconsistent 
with zero bispectrum. 

X^(0) is similar to what several authors have used for quantifying statistical significance of 
non-Gaussianity in the DMR map (Ferreira et al. 1998; Magueijo 2000; Sandvik and Magueijo 
2000). They use only diagonal terms of the covariance matrix; however, the matrix is diagonal 
only on the full sky. As lack of sky coverage correlates one mode to the others, we should include 
off-diagonal terms as well. We did so in equation (12). 
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Fig. 5. — Constraint on the non-linear coupling parameter, /nL) which characterizes non-linearity 
in inflation (Eq.(lO)). The dashed, solid, and dotted lines represent the three different Galactic 
cuts as quoted in the figure. The thick vertical lines plot the measured values of /nl from the 
COBE DMR maps, while the histograms plot those drawn from the Monte-Carlo simulations for 
each cut. 68% confidence limits on /nl are |/nl| < 1-6 x 10'^, 1.5 x 10^, and 1.7 x 10^ for the three 
cuts, respectively. 
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Fig. 6. — Constraint on amplitude of the interstellar dust bispectruni, /dust- The meaning of the 
lines is the same as in figure 5. 
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Fig. 7. — Constraint on amplitude of the interstellar synchrotron bispectrum, /sync- The meaning 
of the lines is the same as in figure 5. 
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Figure 8 compares Xy)Mr(.^) with Xmc(O) ^^^ ^^^ different Galactic cuts. The measured val- 
ues are Xdmr(O) ~ 475.6, 462.7, and 464.8 for the corresponding Galactic cuts, respectively, 
while (xmc(^)) ~ ^^^- ^^ ^^'^ *^^ probability of Xmc(^) being larger than Xdmr(O) *o be 
P (Xmc > Xdmr) = 36.9%, 47.0%, and 49.7%, respectively. 

We conclude that the DMR map is comfortably consistent with zero normalized bispectrum. 
We explain the claimed detection (Ferreira et al. 1998) by a statistical fluctuation as an alternative 
to the "eclipse effect" proposition made in Banday et al. (2000). 

There is no evidence that the scatter of the normalized bispectrum is too small to be consistent 
with Gaussian, in contrast to the claim of Magueijo (2000) based on x^(0) derived from 8 modes. 
To clarify, our analysis does not reject the possibility that the CMB sky is non-Gaussian for only a 
small number of modes; however, in the absence of a theoretical motivation for limiting the analysis 
to a specific set of modes, we choose to treat all the bispectrum modes on an equal footing. Sandvik 
and Magueijo (2000) claim that the non-Gaussianity found by Magueijo (2000) does not spread to 
other modes. This is consistent with our result. 

Incidentally, we plot in the figure 8 the x^ distribution for 466 degrees of freedom, xlee' ^^ 
filled circles; we find that the distribution of Xmc(O) ^^ ^^^y similar to the xlee distribution for a 
smaller cut as expected, while it becomes slightly broader for a larger cut for which the distribution 
of the normalized bispectrum deviates from Gaussian appreciably. Yet, we find that the 20° — 25° 
cuts reasonably retain the Gaussianity of the distribution of the normalized bispectrum. 



5. DISCUSSION AND CONCLUSIONS 

In this paper, we have measured all independent configurations of the angular bispectrum on 
the COBE DMR map, down to the DMR beam size. Using the most sensitive sky map to CMB, 
which combines the maps at 53 and 90 GHz, we test the Gaussianity of the DMR map. 

We find that the normalized bispectrum, Bi^i^iJ {Ci^Ci^Ci,^) ' , gives more robust test of 
Gaussianity than the bare bispectrum, Bi^i^i^. We compare the measured data with the simulated 
realizations, finding the DMR map comfortably consistent with Gaussian. We explain the reported 
detection of the normalized bispectrum at li = I2 = l^ = 16 (Ferreira et al. 1998) by a statistical 
fluctuation. While it is still conceivable that the eclipse effect of the Earth against the COBE 
satellite generates some of the bispectrum (Banday et al. 2000), the DMR data cannot distinguish 
it from the statistical fluctuations. 

We fit the predicted bispectra to the data, constraining the parameters in the predictions, 
which include the primary bispectrum from inflation and the foreground bispectra from interstellar 
dust and synchrotron emissions. We flnd that neither dust nor synchrotron emissions contribute 
to the bispectrum significantly. 

We have obtained a weak constraint on the non- linear coupling parameter, /nL; that charac- 
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Fig. 8. — Testing hypothesis of the normahzed bispectrum, Bi-^i^i.J {Ci-^Cifii.^) ' , being zero in 
the C05-EDMR four-year 53 + 90 GHz sky map. The dashed, solid, and dotted lines represent the 
20° cut, the extended cut, and the 25° cut, respectively. The thick vertical lines plot the measured 
X^(0), while the histograms plot those drawn from the Monte-Carlo simulations. The filled circles 
plot the x^ distribution for 466 degrees of freedom. 
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terizes non-linearity in inflation. We interpret the constraint in terms of a single field inflation as 
follows. According to the analysis of non-linear perturbations on super horizon scales (Salopek and 
Bond 1990), we can explicitly calculate /nl as 

5 f d^lnH \ 

where H is the Hubble parameter during inflation. When applying the slow-roll conditions to an 
inflaton potential y(0), we have (9 In if /9(?i) ~ {d\nV/d4>)/2; thus, /nl is on the order of curvature of 
a slow-roll potential, implying that |/nl| should be smaller than 1 in slow-roll inflation. Therefore, 
the obtained constraint, [/nlI < 1-5 x 10^, seems too weak to be interesting; however, any deviation 
from slow-roll could yield larger |/nl|i bigger non-Gaussianity. 

The next generation satellite experiments, the Microwave Anisotropy Probe (MAP) and Planck, 
should be able to put more stringent constraints on /nl- Komatsu and Spergel (2001) have shown 
that MAP and Planck should be sensitive down to |/nl| ~ 20 and 5, respectively. We find that 
the actual constraint from COBE (figure 5) is much worse than the estimate. This is partly due to 
different cosmology used for the model, but mainly due to incomplete sky coverage; the statistical 
power of the bispectrum at low multipoles is significantly weakened by the Galactic cut. Since 
MAP and Planck probe much smaller angular scales, and their better angular resolution makes an 
extent of the Galactic cut smaller, the degradation of sensitivity should be minimal. Moreover, the 
improved frequency coverage of future experiments will aid in extracting more usable CMB pixels 
from the data. At this level of sensitivity, any deviation from slow-roll could give an interesting 
amount of the bispectrum, and MAP and Planck will put severe constraints on any substantial 
deviation from slow-roll. 

While we have explored adiabatic generation of the bispectrum only, isocurvature perturbations 
from inflation also generate non-Gaussianity (Linde and Mukhanov 1997; Peebles 1997; Bucher and 
Zhu 1997). They are in general more non-Gaussian than the adiabatic perturbations; it is worth 
constraining those models by the same strategy as we have done in this paper. 

We would like to thank Charles L. Bennett, Gary Hinshaw, Misao Sasaki, Licia Verde, and Ed- 
ward L. Wright for helpful discussions. We would like to thank Uros Seljak and Matias Zaldarriaga 
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Japan Society for the Promotion of Sciences. D. N. S. and B. D. W. are partially supported by the 
MAP/MIDEX program. 



A. ANGULAR POWER SPECTRUM AND BISPECTRUM ON THE 

INCOMPLETE SKY 

Incomplete sky coverage destroys orthonormality of the spherical harmonics on the sky. The 
degree to which orthonormality is broken is often characterized by the coupling integral (Peebles 
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1980), 

Ww^^> = I cfh Win)Yil (A) Yir,^. (n) = f Su Y^^^ (n) Yv^> (ii) , (Al) 

where VF(n) is zero in a cut region otherwise 1, and ilobs denotes a soHd angle of the observed sky. 
When Wii/^m' ¥" ^ii'^mm'-, the measured harmonic transform of the temperature anisotropy field, 
aim-, becomes a biased estimator of the true harmonic transform, a[Jj^'^, through 

oo V 

aim = 2^ 2^ OVrn'^ll'mm'- (A2) 

l'=Om'=-l' 

Hence, we must correct our estimators of the power spectrum and the bispectrum for the bias 
arising from incomplete sky coverage. 

First, we derive a relationship between the angular power spectrum on the incomplete sky and 
that on the full sky. Taking the ensemble average of the estimator of the power spectrum, the 
pseudo-Q (Wandelt et al. 1998, 2000), Q = {21 + 1)-^ J2m Wimf, we have 



I' mm' 

« ^CT' E E / ^'^ W{n)Yi*m (n) Yi,m' {n) f d^m W{rh)Yim (m) Y^^, (m) 

m I'm' •' •' 

^^true ^ j ^2^ Wih)YCm (n) J <fih W {rh)Yim (ih) 5^^^ (u - rh) 
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m 

cr^" j'^w{u)Pi{i) 



-Ytruc ^ ^obs 



= Cr^^. (A3) 

In the second equality, we have taken C\f^^ out of the summation over V , as \Wii'mm'\ peaks very 
sharply at / = /', and Cff^'^ varies much more slowly than {Wwmm'l in ^'- This approximation is 
good for nearly full sky coverage. In the third equality, we have used J2i'm' ^I'm' (n) Y^*^, (rh) = 
j(2) ^^ _ j^^_ jj^ ii^Q forth equality, we have used ^^ Yj^ (n) Yim (m) = ^^P/(n • rh). The result 
indicates that the bias amounts approximately to a fraction of the sky covered by observations. 

Next, we derive a relationship between the angular bispectrum on the incomplete sky and that 
on the full sky. We begin with 



{aiimiahmiai-jm-i) - E \4[m;4>^'^|m^/^«i/imimi^/2«>2m^^«3«3m3"^3- (^^) 

all I'm' 

Rotational and parity invariance of the bispectrum implies the bispectrum given by 

(Oiimia/ama^ma) = ^hhh / ^^^ ^hmi («) Yi*^^^ (h) Yil^,^ (h) , (A5) 
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where bi^^i^i^ is an arbitrary real symmetric function, which is related to the angular averaged 
bispectrum, Bi^i^i^. When hY^\ varies much more slowly than the coupling integral, we obtain 

(a/imiazamaaiams) = Y^ ^X^i^v.^ Y^ j rf^n Y^l^^,^ (A) Y^r^, (A) F^f^, (A) 

all V all m' '' 

X j <fhr W{n,)Yv^ra', i.i^i)yU, (Ai) 
X j <fu2 W{u2)Yi,^ra', ("2) >i:„, (na) 

^ b\llxld'nW{h)Yil„,^{h)Yi;^Jn)Yi;^^{h). (A6) 

Then, we calculate the angular averaged bispectrum, Bi-^i^i^ (Eq.(l)). By convolving equation (A6) 
with the Wigner-3j symbol and using the identity (3), we obtain 

-1 



(Bi^mJ « b 



, _ true / ^ f h h h \ 

nM',J - '"'I'^Zay (2li + l)(2^2 + l)(2^3 + l) \ 0^ 
X ^ / (fm Yi^rni (m) yi2„2 (m) Yi^ms (m) 

all m 

X I d'n W{n)Yil„,^ (n) y,*„^ (A) y,*„3 (A) 

-1 



6; 



true / (2ll + l)(2/2 + l)(2^3 + l) / /l h h \ 

''I'^'sV 4^ 1^ J 

/d YJi f d n 
~4^ / 17 ^(")^'i ("^ ■ ") ^'2 (m ■ n) Pi3 (m • n) 

.true /(2/l + l)(2/2 + l)(2/3 + l) /^ ^1 ^2 ^3 ^ f^obs 



4tt \ / 47r 

where we have used the identity. 



itrue ^^obs 



'^^ P,,(x)P,,(x)P,3(x) = f p' p' p' ") . (A8) 

Thus, the bias for the angular bispectrum on the incomplete sky is also approximately given by a 
fraction of the sky covered by observations. 
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